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1. Introduction 

Dissipative processes in isotropic and homogeneous Friedmann-Lemaitre-Robertson- 
Walker (FLRW) universes are restricted to scalar dissipation by the spacetime 
symmetries. Such scalar dissipation may be analyzed via the relativistic theory of bulk 
viscosity (see [|fl |2| and references cited there). Relativistic thermodynamics considers 
processes which remain close to equilibrium, so that the transport equation is linear 
in the bulk viscous stress. The original and often used theory of Eckart, and a similar 
theory due to Landau and Lifshitz, are non-causal (they admit superluminal signals) 
and their equilibrium states are unstable ||. The theory of Israel and Stewart ||, ||], 
and related theories ||, overcome these pathologies, since they are causal and stable 
under a wide range of conditions J7|. 

There are dissipative processes in cosmology which are not close to equilibrium. For 
example, inflation driven by a viscous fluid necessarily involves a bulk viscous stress 
that exceeds the equilibrium pressure |], |J. The reheating era at the end of inflation 
also involves far-from-equilibrium dynamics. In ||, reheating was analyzed via causal 
thermodynamics. In that model, the bulk viscous stress is small, and the dominant non- 
equilibrium effects arise from particle creation. It is also possible to consider models 
with an increased contribution from the bulk viscous stress, beyond the linear regime. 

In order to treat these and other dissipative processes that do not remain close to 
equilibrium, one requires a nonlinear generalization of the standard theories. Recently 
a nonlinear theory was developed in M, via a simple phenomenological generalization 
of the Israel-Stewart theory. The nonlinear theory reduces to the causal theory in the 
linear regime, and the second law of thermodynamics is built into the theory. There is 
consequently an upper limit to the bulk viscous stress, as in nonlinear generalizations 
of the non-relativistic heat flow equation || and shear viscous stress equation ]10[ . 

One can define at each event a local reference equilibrium state with energy density 
p, pressure p, number density n, specific entropy s and temperature T. Then the energy- 
momentum tensor of the fluid for the case of scalar dissipation (i.e., no energy or particle 
flux and no anisotropic stress) is 

T a p = (p + V + TVjUaUp + (p + n)# a/3 , 

where IT (< 0) is the bulk viscous stress, u a is the four- velocity with respect to which 
p, p and n coincide with the local equilibrium values (in FLRW spacetime, this will be 
the preferred four- velocity), and g a/ 3 is the metric (specialized to FLRW below). The 
effective non-equilibrium specific entropy is |6|. ^ 

where r(p, n) is the characteristic timescale for linear relaxational effects (the crucial 
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thermodynamic parameter in Israel-Stewart theory that ensures causality), and £(p, n) is 
the linear bulk viscosity. The evolution of the equilibrium and effective specific entropies 
is given by || 
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where H = ^V a u a is the Hubble rate and is a characteristic timescale for nonlinear 
effects. The second law of thermodynamics imposes an upper limit on the bulk stress 



mi < 



c 



(4) 



(Note that in the linear theory, there is no upper limit on |II| that is built into the 
theory - |H|/_p <C 1 is implicitly imposed a priori since the theory is linear.) 
The nonlinear transport equation for II is M 




11(1 + 3nH) 
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Israel-Stewart theory (in its full, 'non-truncated' form |7|, [I|]) is recovered when r, 
If in addition r = 0, then the transport equation reduces to that of Eckart. 
In a flat FLRW universe, Einstein's equations aret 

H 2 =|p, (6) 
H + H 2 = - |(p + 3p + 3II), (7) 

so that by (5) and (7), we get a second-order evolution equation for H when r ^ 0. 
In order to find this evolution equation for nonlinear dissipative processes in a flat 
FLRW universe, we need to specify the thermodynamic parameters (, r and r*, and the 
equations of state relating p, p, T and n. These will depend on the physical situation 
being analyzed, and all quantities except r* follow from standard thermodynamics. In 
the nonlinear theory of ||, the microscopic interactions governing far-from-equilibrium 
behaviour are unknown and encoded into the phenomenological parameter r*. As a first 
step towards a more complete analysis, we follow || and postulate here a simple ansatz 
for r*, which allows us to investigate the qualitative nature of the effects of nonlinearity: 



— k T , 



(8) 



f We use units with 8ttG = 1 = c and &b = 1. 
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where k is a dimensionless constant. This assumption is a mathematical simplification, 
but will allow us to find some overall features of nonlinear viscosity (whose detailed 
physics is not known). The weakly nonlinear regime is then characterized by A; 2 <C 1, 
and we will linearize the evolution equation for the Hubble rate below. 

The linear relaxation time is related to the bulk viscosity by the physical relation 



11 



(9) 



v 2 (p + p) 



where v is the dissipative contribution to the speed of sound V (< 1), so that 
V 2 = c 2 + v 2 , with c s the adiabatic contribution. This relation follows from an analysis 
of the propagation of small disturbances in a dissipative fluid. A simple form for the 
bulk viscosity is often taken to be ( oc p q l 2 (q constant), so that by (6) we have 

C = aH\ (10) 

where a (> 0) is a constant. Although this is a mathematical ansatz introduced for 
simplicity, it does approximate the physical form of ( for certain fluids \T2\. 



We also need to specify the equations of state. For the pressure, we will take the 
simple linear barotropic law 

p=(7-l)p, (11) 

where 7 is a constant. Typically 7 ~ | (radiation-dominated expansion) or 7 « 1 
(matter-dominated expansion). (Note that the limiting case of dust, when p is exactly 
and 7 is exactly 1, is strictly ruled out, since the temperature is zero for dust and there 
can be no bulk viscous stress.) By (JTl~|), c 2 = 7 — 1 is a constant, and we will assume 
that v is also constant, subject to the causality constraint 

v 2 <2-j. (12) 

This is a mathematical simplification, but it should not be unreasonable if (|TTj) holds. 

By equations (6) and (g)-(|n]), we can use the above thermodynamic assumptions to 
rewrite the upper limit (|) on the bulk stress, and to express r* in terms of the Hubble 
rate: 



ITT. C (V\ 2 



7(2-7) 



A' 2 



P, (13) 



It follows from (14) that in the case q — 1, the nonlinear (and linear) timescales are 
directly determined by the expansion timescale H" 1 . 
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Finally, for the temperature, we will consider two cases: (a) barotropic temperature, 



T = T(p), so that [11 



T oc p^~ l)h , 
and (b) ideal-gas temperature, so that 
p = nT . 



(15) 



(16) 



Note that ( |T5D and (|T^) can only hold simultaneously in equilibrium O, [T3] ]. 

Equations (6)— (|TT]) together with ([TB]) and the number conservation equation 
h + 3Hn = give the evolution equation for case (a) as || 



— IH + 3HH + 




1l 



7 



H A 1 



[2H + 3 7 # s 
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(17) 



In case (b), the ideal-gas behaviour (0) leads to the new and simpler evolution equation 



2k 2 

3r^lp 



HH - 2H< 
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+ ^^H 3 ~ q (2H + 3 7 if 2 ) = 
2a 



Considerable work has been done using the linear non-causal model of bulk viscosity 
T2| . [l| for references) , including the case of cosmic strings [[14], [15|, [16], |17| . Our results 



(see 

represent a generalization to the nonlinear case of previous dynamical analysis in the 
linear causal theory |12, 13, 18, 19, 2C, 21], and also a generalization of some of the 
results in ||, where exact de Sitter solutions for any q, and a power-law solution with 
q = 1 were given for equation (17). Here we find new solutions with q — 1 in the weakly 
nonlinear regime for both (17) and (18). We also investigate asymptotic solutions and 
their stability for q 1. In principle our solutions and their properties could be used 
to model far-from-equilibrium processes, such as viscous fluid inflation and reheating. 
However, we will not pursue such applications here. 
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2. Models with barotropic temperature 

2.1. Solution for q = 1 

When the bulk viscosity index satisfies q = 1, we have from fllCf ) and (14) that ( oc H 
and r, r* oc Assuming that fc 2 < 1 and H/H 2 is bounded, we can isolate H in 

(17) to first order in k 2 : 



HH 



l-2 7 , 2fc 2 \ rV2 , ^ , 7^ 2 , 2 7 A; 2 1 ,, 2 



+ I iF + 3 1 + — + 



H H 



+ 



97/I 7V 2 2 , 7/c 2 ' 



+ 



2 \2 a 
With the change of variables 

H = y n , t = 0/(3m), 

where 



# 4 = 0. 



(19) 



1 (I-7) , 2k 2 7 v 2 2 7 fc 2 
- = 1 , m = H 1 

n 7 a a a 



equation (19) becomesf 



do 2 J 



■2m 



where 



dO n+V 
a (a + 27« 2 — 2at> 2 ) 



+ 



71; (a + 71) 2 - «7V ) 



4 (a + 7t> 2 ) 2 
Then, making the nonlocal transformation |20] 

z = _L- 3/ «+ 1 , d v = y n de, 

n + 1 

equation (|20|) becomes linear: 

G? 2 2: 

— + — + /fc = 0. 
drj z drj 



[a + 7W 2 )" 



(20) 



A: 2 + 0(A; 4 ). (21) 



(22) 



(23) 



From (|2lD it is easy to check that 1 — 4/3 > 0, so we obtain the general solution of (19) 
in parametric form 

H(rj) = [{n + 1) {C + exp (A+tj) + C_ exp (A_r?)}] n/(n+1) , 

o??7 



^ =t0+ (i)/ 



(24) 
(25) 



f Note that n — — 1 is ruled out since fc 2 /a is non-negative. 



1 /2n 

where A± = [— 1±(1 — 4/3) j/2 are real, and t , C± are arbitrary integration constants. 
When either of C± vanishes, we obtain one-parameter families of solutions, which have 
the form: 

2 7 (v 2 + k 2 ) 



H±(t) 



it - 1 ) 

^0 



it - 1 ) 



if a^ao = 
H + = Hq if a 



2v 2 - 1 

- a , 



where 



2 

37, 



m ± {m 2 - [to + (7/a - 2) w 2 ] (1 + 27fc 2 /a)} 



1/2 



to + (7/a — 2) 



^0 



a + 2 7 A; 2 
3a 7 m 



(26) 
(27) 

(28) 
(29) 



and Hq is an arbitrary positive constant giving a de Sitter solution (with v > l/\/2), 
which agrees with the linearization of the de Sitter solution given in ||, equation (41). 

The parameter v_ is smaller than the perfect fluid value: < z/_ < 2/ (37) < 1, and 
for small dissipative contribution to the speed of sound we get v- ~ 2 (3 7 ) _1 
Also vq < (3 7 ) _1 < 1 for 1 < 7 < 2, while —00 < u + < 00. Thus, (27 
singular solutions without particle horizons. 

Evaluating the local-equilibrium entropy s from (2), we find 



1 - y/2v 
) contains 



s(t) = s + 7 



2 W7 fl 3 



n T 



for the solution given by (26), and the same for (27 
effective specific entropy s e s is given by ([[]): 



(t - tof^- 2 ^ , (30) 
but with uq instead of v±. The 



s e s(t) = s + 



(3u 2 ± ) 



2 n!/7^3 



n T 



7 



1 /37z/±-2' 



2v 2 7 V 3i/± 



to 



,3^-2/7 



Hence the requirement s c q > is fulfilled for 

3i/±7 - 2 



3z/j 



< V2vy . 



.(31) 



(32) 



For the solution (27-), the restriction (|32|) and the solution (^Tj) are the same, but u± 
is replaced by z/(j. Provided that vq > 0, then this restriction requires v > l/y/2. 
The second restriction, arising from positivity of the entropy production is 



3^±7 



2 7V 2 
3i/± ~ < l 2 " 



(33) 



for the solution (26), and the same but replacing v± by uq for (27). Thermodynamics 
restrictions (|3"2"D and (B3) are the same one for v = \[2k 2 . 
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We are able to find a third solution. Provided the parameters satisfy the constraint 

7 [m + ( 7 /g - 2) v 2 ] = a 2 (l - 7) + 2^ak 2 

8am 2 ' [a(2- 7 ) +4 7 fc 2 ] 2 ' 1 J 

the general solution of (19) takes an explicit form, and we get for the scale factor 

, > , (n+2)/3mn 

a(t) =a \t-t \ n+1 + K , (35) 



where K is an integration constant. We have checked that the constraint (|34|) is satisfied 
for suitable values of the parameters. 

As examples, we apply the two thermodynamical restrictions ([32]) and (|33|) and 
the causality constraint fljjl) to the solutions (26) and (27), for a radiation-dominated 



expansion, 7 = |. Causality demands in this case < |. 

• For the solution (26) we choose v 2 = | and a = 1. Causality and the two 
thermodynamical restrictions are fulfilled for any k 2 <C 1. By (28), it follows that 
for small k (k < 0.32), we have u + > 1, so that these solutions are always inflationary, 
and i/_ < 1, corresponding always to non- inflationary expansion. The maximum value 
of v + is I (5 + vl3), and the minimum value of z/_ is |(5 — VT3). Thus a ~ t N where 
0.17 < N < 1.08. For small a (perfect fluid limit) we get v + « 2(37) _1 < 1 and 
V- ~ (2/3-y)(k 2 /(k 2 + v 2 )) 1. Thus the solution with z/ + corresponds to the perfect 
fluid limit solution, while the other corresponds to a new solution that does not exist in 
this limit. On the other hand, for large a we get 

2 1 



37 1 =|= \[2v 

thus, v± may be very large for v 2 close to 1/2 (expanding solutions require v 2 < 1/2). 



For the solution (27-), we choose v 2 — |, and by (29) 

^0 = ^(l-|A: 2 ) +0(A; 4 ), a = § (2 + 3k 



The solution is not inflationary, and causality and thermodynamic requirements are 
always fulfilled. 

The requirement that H/H 2 remains bounded means that not all the two-parameter 
families of solutions (25) of equation (19) are also solutions of (17), but it is a 
consequence of the thermodynamic restriction given by positivity of the entropy 
production. So, from (|J), (6)-([5|) and (|TT|), we find that H/H 2 is bounded: 

H 2 ~2k 2[V k) - 
Solutions that satisfy this behave as follows ]H| : 

The evolution begins at a singularity with a Friedmann leading behaviour of the form 
(26-) for a 7^ «o or (27-) for a = «o- Then the expansion becomes 
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1. asymptotically Friedmann, like (26+), for a < a ; 

2. asymptotically de Sitter, like (27+), for a = a ; 

3. divergent at finite time, with leading behaviour (26+), for a > «o- 



2.2. Asymptotic solutions for q > 1 and stability 



For q > 1 it is easy to check that (17) admits a solution whose leading term is 2/(3jt) 
when t — ■> oo. To study its stability, we make the change of variables 

H- 1 = t[u(^)] 7/(7 - 1} , t q ~ l 



when 7 > 1. Then (17) takes the form 



u"4 
+ 



671/ 



-u 
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3 (7 - 1) v 2 r 



7 (g-2)/( 7 -l) + 



(q-lyaR 



z, 
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z 
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z \ q — 1 



u (79-7-l)/(7-l) _ ^Z M (79-2 7 -l)/(7-l) 



7-1 

(ff-l) S 
= 0, 



9 2fe 2 -^ 2 + V + i)/a-7) + 3^±^ M i/d-7) 



i? 



where a prime denotes d/cfe and 

' 2fc 2 ' 
,37f 2 , 

We may rewrite (36) in the form 



k 2 

R=l--z 

v 



u 



7/(7-1) 



1 + 



7# 



7 ((7 — 1) zw' 
7 — 1 u 



u 



r 



d \i„./2 
dz 



u' 2 + U(u,z) =D(u,u',z) 



(36) 



(37) 



(3f 



When z — > 00, equation (36) has a constant solution w = (§7) 7 7 , which is the 
unique minimum of U(u, z). Thus, considering that u lies in a neighbourhood of uq, we 
get from ([J7]) that 



R^l- 



7 



1/7 



v 



zu 



7 — 1 \37/ 

We can see how R — > 1 when z — > 00 through a linearization of (36): 

(u - u ) 



(39) 



u" + 2(g - 1)£V + n- 







(40) 



where Q = 3(j /a)v 2 (3j/2) q 2 (q — 1) 2 . By the change of variable u = uq + u> exp[(l 
?)Oz], equation (|40|) becomes 

ft 



w" + 



-(g-l) 2 ft 2 + 



io = 0, 



(41) 
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which can be solved in terms of confluent hypergeometric functions. It has the following 
leading behaviour for z — > oo: 

w{z) ~ C lZ 1/2{1 - q) exp [{q - l)Qz] + C 2 z 1/2{q - 1) exp [(1 - q)Qz] , (42) 

where C\p, are arbitrary integration constants. Thus zv! ~ z 1 ! 2 ^^ — > for z — > oo and 
q > 1, and we get the expansions 



t/"(u,«) = 
for 7g 7^ 2 or 2 + 7, 

C/(u,z) = 

for 7g = 2, 

C/(u,z) = 
for 7g = 2 + 7, and 

D(u, u' z 



3( 7 -l)V 
(q - if a 

3( 7 - l)v 2 
(q-l) 2 a 

3 (7 - 1) t; 2 
(q- if a 

67V 2 



(79-2)7(7-1) 



37M 



(79-2- 7 )/(7-l) 



iq - 2 2 ( 79 - 2 - 7) 





vz 2 ) 







M 2 



1^7/(7-1) _ ^2 ^ 

7 2 m 



M 7( 9 -2)/(7- V2 + O Q . 



,(43) 



(44) 



(45) 



(46) 



a(q - 1)_ 

For mathematical completeness, we will consider also the case 7 = 1. We make the 
change of variables 

H' 1 = texp [-u(z)] , z = t q ~ l , 

in (17), which then takes the form 



u 



1 (3f_ 
2{Ra 

3v 2 



e -(9-l)« + 



(q-2)e- u + 3(^ + 1 



u 



z J q — 1 



1 



4i 



2e~ u + 9 



fc 2 -t; 2 1 
R + 2 



e "" 12 F = °- (47) 



When 2; — > 00, equation (47) has a constant solution w — hi|. Again, we may write 
(47) in the form (|38D , and find that uq is the unique minimum of U(u, z). Thus we get 



for q 7^ 2, 



C/(u, z) 



C/(u, z) 



3v 2 



4a (g- 1)' 



q-1 



-{q-l)u _ 



,-(«-2)« 



q-2 



3^ 
4a 



( 2e_U - 3m ) \ + 



T + Ofi) > (48) 



(49) 
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for q = 2, and 



D(u, u, z) 



3v 2 



(q-l)u 12 




) 



(50) 



2a (q - 1) 



e 



As U(u,z) has a unique minimum at u for any g ^ 1, and D(u,u',z) is negative 
definite for q > 1, we find that solutions with leading Friedmann behaviour a ~ t 2 /( 3 T) 
when t — > oo are asymptotically stable for q > 1. 

3. Models with ideal-gas temperature 

3.1. The de Sitter solution 

Using the ideal gas equation p = nT instead of the power-law for the temperature one 
obtains (18). Setting H = Ho in (18) then gives 



The upper bound (13) on the bulk stress implies the limit k < v. Thus de Sitter solutions 
may exist. On the other hand the specific entropy may be calculated and one obtains 



Linear analysis around the fixed point (0, Hq) gives the following information: for 
1 — ux < q < 1, it is an asymptotically stable node; for q > 1, it is a saddle point; 
for q < 1 — u>i, it is an asymptotically stable focus, where 



For q = 1, de Sitter inflation occurs for H arbitrary, provided that a = v 2r y/(v 2 — k 2 ). 
3.2. Solution for q — 1 

When q — 1, power-law solutions a ~ t exist if, by (18), 



(for 79 ^1) 






(51) 



For positive entropy generation we require 




(52) 
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but it is easy to show that this restriction is always fulfilled for ( |51~D with arbitrary v, k 
and a. From ([!]) and (2) with 7 7^ 1, we find 

(3 7 iV - 2) 2 



s(t) = s + 



37^-2 

. 7-1 . 



hit, s eff (t) = s(t) 



18N 2 v 2 ^ - 1) ' 



so that positive s requires N > 2/3 r y - which is always fulfiled for expanding solutions, 
as one can see from (|5l"l) . These solutions (N > 0) require 



k 2 > 



a — 7 



In order to get an inflationary power-law solution, it is necessary that a > a m \ n where 

w 2 7(3 7 - 2) 



o:, 



3^(v 2 - k 2 ) + 2k 2 



For q — 1, and assuming that k 2 /v 2 = e 2 C 1 and H/H 2 is bounded, we can isolate 
E in (18) to first order in e 2 : 



EE 

9 2 
+ - W7 



37V 2 



O: 



1 + 2e 2 )E 2 E + 2 



1 



a 



l + ^(l + e 2 ) 



E 4 = 0. 



(53) 



With the change of variables E = y n , t = a^/[3 7 f 2 (l + 2e 2 )], where n = a/ '(2k 2 — a) 
equation (53) takes exactly the same form as (PUD, with 

e 2 (7 - a + 7e 2 ) 



(54) 



7(l + 2e 2 r 

Then, making the same nonlocal transformation ( p2|) as before, equation (^) is again 
brought to the linear form fl2~3"|). From ([54]) we find that 1 — 4(3 > 0, so we obtain the 
general solution of (53) in parametric form 



H{v) = {{n + 1) [C+ exp (A+r?) + C_ exp (A_?y)]} 



t + 



a 



3-fv 2 (1 + 2e 2 



drj 



(55) 
(56) 



1 /2i 

where A± = [—1 ± (1 — 4/5) ]/2 are real, and to an d C± are arbitrary integration 
constants. When either of C± vanishes, we obtain one-parameter families of solutions, 
which have the form: 



H±(t) 



it - 1 ) 



if a ^ a = 7 1 + —z 



(57) 
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where 




^=Urh-r^- ( 6 °) 



The specific entropy s is given by fl3U| ) and s e g by (|5T|) , and the restrictions (^) and (R3) 
also apply here, with v± and given by (59) and (60), and satisfying similar restrictions 
to the barotropic temperature case. Note that «o is the value of a which arises from 
the requirement that the solution of (53) is de Sitter, and it can be also obtained as the 
first order expansion of the expression ckq = 7/ (1 — e 2 ) that arises from (18). 
On the other hand, provided the parameters satisfy the constraint 

1 - 8^e 2 + 16?/; 2 e 4 . . 

(61) 



1 + 2ip + e 2 - 4?/> 2 e 2 



where <fi = 7/a and ip = v 2 /a, the general solution of (53) takes an explicit form, and 
we get for the scale factor 

a(t) = a \\t-t \ n+1 +K\" , (62) 

where K is arbitrary and 

4k 2 - a . , 

^VrTrf)' (63) 

We have checked that the constraint ( |5Tj ) is satisfied for suitable values of the parameters. 
For instance, if e < 0.1, then <p < 1. This constraint also reduces, in the limit e — > 0, to 
the constraint on the parameters found in a previous paper for ideal-gas temperature 
in the linear theory [^TJ. It must be remarked however, that the form of the solutions 



is quite different. This shows that, in general, the solutions are not analytic in k for 
k — > 0. In other words, this limit is singular. 

The requirement that H/H 2 remains bounded means that not all the two-parameter 
families of solutions (56) of equation (53) are also (approximate) solutions of (18). Those 
that satisfy it behave as follows: The evolution begins at a singularity, with Friedmann 
leading behaviour as in (57-) for a 7^ «o, or (58-) for a = «o- Then the expansion 
becomes 

1. asymptotically Friedmann, as in (57+), for a < «o; 

2. asymptotically de Sitter, as in (58+), for a = a ; 

3. divergent at finite time, with leading behaviour (57+) for a > ato- The 
relationships for s and s e s and the the thermodynamical constraints ( |32|) and ( |33"D in 
the limit t — > 00, are the same, but with v± replaced by fi(n + 1). 
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3.3. Asymptotic analysis for q > 1 and stability 

For q > 1 it is easy to check that (18) admits a solution whose leading term is 2/(3jt). 
To study its stability, we make the change of variables 

H- 1 = tu{z) , t q - x = z , 

in (18), which then becomes 

/6V , q\ u' 



u" + —-u*-' + - 
V an z 



37 ^^-'-^"- 2 ) + ^^-0, (64) 



{q-l)aRz\ 2 J 2(q-l)z 2 u 



where 



R = 1 - e 2 + — [u + (q - 1) zu'] 



37 

We rewrite (64) in the same form (]38f) as before. When z — > oo, equation (64) has a 
constant solution uq = |7, which is the unique minimum of U(u,z). Thus, considering 
that u lies in a neighbourhood of uq, we get 

R « 1 + — e 2 (g - 1) zu' . 

37 

We can see how i? — > 1 when z — ► oo through a linearization of (64), which again 
leads to (f40|), with Q = 3(^/a)v 2 (q — l)~ 2 (37/2)' 7 ~ 2 . By the same change of variable 
it = Uq + wexp[(l — equation (|40| ) is again brought to the form (^TJ). As before, 

this can be solved in terms of confluent hypergeometric functions, and has the same 
leading behaviour (fi2|) for z — > oo. Thus zu' ~ ^V 2 ^-?) — > for z — > oo and g > 1, and 
we get the expansions 



U{u, z) =- -g- 



(q- 1) a 



? 2(9-1 



z Vz 2 



As C/(u, z) has a unique minimum at u for any q ^ 1, and D(u, u', z) is negative definite 
for g > 1, we conclude that solutions with leading Friedmann behaviour a ~ 

t 2/(3 7 ) when 

t — > oo are asymptotically stable for q > 1. 
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4. Concluding remarks 



With the aim of treating dissipative processes that do not remain close to equilibrium, 
we have used the nonlinear generalization of the standard causal theories proposed in 
||. Barotropic and ideal-gas forms for the temperature have been used to close the 
dynamical equations of the universe. 

The evolution of the universe is qualitatively similar in both cases, and we can 
summarize the main results as follows. Singular solutions are found for q — 1 with 
a power-law and inflationary behaviour in the limit t — > oo. The initial and final 
behaviour of the solutions follows closely that of the linear theory with barotropic 
temperature, in the 'full' causal as well as the 'truncated' causal and the noncausal cases. 
The logarithmic behaviour found for ideal-gas temperature in no longer appears 
when the nonlinear transport equation is considered. Thermodynamic restrictions, 
arising from positivity of the entropy, have been also derived. The evolution begins 
at a singularity with a Friedmann leading behaviour, and the future dynamic is 
asymptotically Friedmannian, asymptotically de Sitter or divergent at finite time, for 
specific values of the bulk viscosity constant a. 

The behaviour of particular de Sitter solutions has already been analyzed in || , for 
barotropic temperature, and it is also derived here for the ideal-gas temperature. For 
q < 1, a stable exponential inflationary phase occurs in the far future when 1 < 7 < 2. 
The condition a = ao is required if q — 1, and such behaviour is unstable for q > 1. 
Our results show that occurrence of viscosity-driven exponential inflationary behaviour 
depends mainly on the equation of state, rather than on the thermodynamical theory 
employed. 

For q > 1 the perfect fluid behaviour a ~ £ 2 /( 3 t) when t — > 00 is asymptotically 
stable, because the viscous pressure decays faster than the thermodynamical pressure. 
However, if q — 1 and a < a , both pressures decay asymptotically as t~ 2 , and the 
exponent becomes v + . The perfect fluid behaviour becomes unstable if q < 1. 

In essence we have analyzed the existence of specific behaviour for a dissipative 
universe not close to equilibrium in terms of the thermodynamic parameters related to 
r, q and the equations of state. 

A final comment concerns the robustness of our results relative to the spatial 
geometry and the source of the gravitational field. Firstly, the question arises as to 
what the effect of non-zero spatial curvature in the FLRW metric will be. In the 
case of non-causal bulk viscosity, where II is algebraically determined by H, it is 
relatively straightforward to investigate this question [TJ], [15|, [L6|, [T7|]. By contrast, in the 
causal theory (linear as well as nonlinear), II is not algebraically determined by H but 
satisfies an evolution equation that couples it differentially to the expansion. Curvature 
introduces the scale factor a = exp / Hdt explicitly into the Friedmann equation, and 
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makes it more difficult to decouple the equations. Thus the effect of curvature in the 
causal case will be far more difficult to determine in general. (In flTB |, curvature is 
considered in the linear causal theory, but with highly simplified equations of state.) 
This is a subject for further investigation. 

Secondly, it could be interesting, but in general very difficult, to determine the effect 
of introducing another fluid, without bulk viscous stress. In the non-causal theory, this 
question is much less difficult to investigate JT5], |TB[ . In the causal theory, this is another 
subject for further investigation. 



Acknowledgments 

VM thanks the Spanish Ministry of Education for partial support under grant PB94- 
0718. LPC and ASJ thank the University of Buenos Aires for partial support under 
project EX-260. 



References 



Maartens R 1995 Class. Quantum Grav. 12 1455 

Zimdahl W 1996 Phys. Rev. D 53 5483 

Lindblom L and Hiscock W A 1983 Ap. J. 267 384 

Israel W 1976 Ann. Phys. 100 310 

Israel W and Stewart J M 1979 Ann. Phys. 118 341 

Jou D, Casas- Vazquez J and Lebon G 1996 Extended Irreversible Thermodynamics (Berlin: 
Springer) 

Hiscock W A and Lindblom L 1983 Ann. Phys. 151 466 
Maartens R and Mendez V 1997 Phys. Rev. D 55 1937 
Zimdahl W, Pavon D and Maartens R 1997 Phys. Rev. D 55 4681 
Bidar H, Jou D and Criado-Sancho M 1996 Physica A 233 163 

Maartens R 1996 in H Rund Workshop on Relativity and Thermodynamics ed Maharaj S D (Natal 



University, to appear) (astro-ph/9609119) 



Pavon D, Bafaluy J and Jou D 1991 Class. Quantum Grav. 8 347 
Mendez V and Triginer J 1996 J. Math. Phys. 37 2906 
Barrow J D 1986 Phys. Lett. 180 235 
Barrow J D 1988 Nucl. Phys. B 310 743 

Barrow J D 1990 in The Formation and Evolution of Cosmic Strings ed Gibbons G W, Hawking 

S W and Vachaspati T (Cambridge: Cambridge University Press) 
Turok N 1988 Phys. Rev. Lett. 60 549 

Coley A A, van den Hoogen R J and Maartens R 1996 Phys. Rev. D 54 1393 
Chimento L P and Jakubi A S 1997 Class. Quantum Grav. 14 1811 
Chimento L P 1997 J. Math. Phys. 38 2565 

Chimento L P, Jakubi A S and Mendez V 1997 Int. J. Mod. Phys. D in press 



